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ABSTRACT 



We revisit the possibility and detectability of a stochastic gravitational wave 
(GW) background produced by a cosmological population of newborn neutron 
stars (NSs) with r-mode instabilities. The NS formation rate is derived from both 
observational and simulated cosmic star formation rates (CSFRs). We show that 
the resultant GW background is insensitive to the choice of CSFR models, but 
depends strongly on the evolving behavior of CSFR at low redshifts. Nonlinear 
effects such as differential rotation, suggested to be an unavoidable feature which 
greatly influences the saturation amplitude of r-mode, are considered to account 
for GW emission from individual sources. Our results show that the dimension- 
less energy density Qqw could have a peak amplitude of ^ (1 — 3.5) x lO"*^ in 
the frequency range (200 — 1000) Hz, if the smallest amount of differential rota- 
tion corresponding to a saturation amplitude of order unity is assumed. However, 
such a high mode amplitude is unrealistic as it is known that the maximum value 
is much smaller and at most 10^^. A realistic estimate of flaw should be at least 
4 orders of magnitude lower (~ 10^^^), which leads to a pessimistic outlook for 
the detection of r-mode background. We consider different pairs of terrestrial 
interferometers (IFOs) and compare two approaches to combine multiple IFOs 
in order to evaluate the detectability of this GW background. Constraints on 
the total emitted GW energy associated with this mechanism to produce a de- 
tectable stochastic background (a SNR of 2.56 with 3-year cross correlation) are 
~ 10"^ MqC^ for two co-located advanced LIGO detectors, and 2 x 10"^ MqC^ for 
two Einstein Telescopes. These constraints may also be applicable to alternative 
GW emission mechanisms related to oscillations or instabilities in NSs depending 
on the frequency band where most GWs are emitted. 
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Subject headings: gravitational waves - stars: neutron - stars: formation - super- 
novae: general 



1. Introduction 

A stochastic gravitational wave background (SGWB) is a target for gravitational wave 
(GW) interferometers (IFOs). It could have two very different origins. It may result from a 
large variety of cosmological processes developed in the very early universe, such as ampli- 
fication of quantum vacuum fluctuations, phase transitions, cosmic strings, etc. (see, e.g., 
Maggiore 2000; Buonanno 2003 for reviews). Additionally, an astrophysical GW background 
(AGWB) is expected to be produced by the superposition of a large number of unresolved 
sources since the beginning of star formation (Schneider et al. 2000; Regimbau & Mandic 
2008). There have been a host of literatures dedicated to the studies of various AGWB 



sources, such as core collapse supernovae (CCSNe) (IBlair fc Ju 19961 ) . leading to the forma- 
tion of neutron stars (NSs) (Coward et al. 2001; Howell et al. 2004; Buonanno et al. 2005) 
or black holes (BHs) (Ferrari et al. 1999a; de Araujo et al. 2004; Pereira & Miranda 2009), 
phase transitions in NSs (Sigl 2006; de Araujo & Marranghello 2009), coalescing compact 
binaries consisting of NSs and/or BHs (Schneider et al. 2001; Farmer & Phinney 2002; 
Regimbau & de Freitas Pacheco 2006a; Regimbau & Chauvineau 2007), magnetars (Regim- 
bau & de Freitas Pacheco 2006b) and population HI stars (Sandick et al. 2006; Suwa et al. 
2007; Marassi et al. 2009) among others. 

In this paper we revisit the possibility that the r-mode instabilities in newly born NSs 
could form a SGWB. NSs, having long been considered to be likely observational sources 
for GW detection, emit gravitational radiation in a number of ways, for example, through 
CCSNe, inspiralling com pact binaries, rotating deformed stars, oscillations and instabilities 



( lAndersson et al. 20101 ). First postulated more than ten years ago (Andersson 1998; Fried- 
man & Morsink 1998), the r-mode instability has been attracting increased attention due 
to the fact that it is driven unstable by GW emission and it can be active for a wide range 
of core temperatures and angular velocities (Lindblom et al. 1998; Andersson & Kokkotas 
2001; Andersson et al. 2003). Early estimate indicated an energy equivalent to roughly 1% 
of a solar mass is radiated in GWs (Andersson et al. 1999) as an initially rapidly rotating 
star spins down. This led to an expectation of a SGWB produced by a cosmological popu- 
lation of young rapidly rotating NSs with closure density h'^flgw peaking at ~ 10~^ of t he 



present-day critical energy density of the universe (lOwen et al 19981 : iFerrari et al. 1999bl ) . 
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The most important aspect of the r-mode intabihty is the largest amphtude (often 
called the saturation amplitude a) that the perturbation can grow to. This maximum 
amplitude determines how fast the NS spins down and whether the associated GW emission 
will be detectable (either in terms of single event or a stochastic background). In Owen 
et al. (1998) and then Ferrari et al. (1999b) it was taken to be of order unity (there 
were no estimates of this maximum at that time). Later more detailed (both analytical 
and numerical) studies seriously questioned the potential of the instability as well as the 
efficiency of GW emission (Rezzolla et al. 2000; Ho & Lai 2000; Lindblom et al. 2000; 
RezzoUa et al. 2001a, 2001b; Lindblom & Owen 2002). On the one hand, it was suggested 
that energy transfer t o other stellar inert i al modes can s i gnificantly reduce the saturation 



amplitude of r-mode (ISchenk et al. 20021 : iMorsink 20021 : iBrink et al. 20041 ). Arras et al 



(2003) tested nonlinear coupling between stellar inertial modes and revealed much lower 
values of saturation amplitude {a ~ 10~^ — 10~^). Then a specific resonant three-mode 
coupling between the r-mode and the pair of fiuid modes was identified to be responsible 
for the catastrophic decay of large-amplitude r-modes and a perturbative analysis of the 
decay rate suggest ed a maximum dimensionless saturation amplitude ttmax < 10~^ — 10^^ 



(ILin fc Suen 20061 ). More recently Bondarescu et al. (2009) examined the 3-mode coupling 
between the r-mode and two other inertial modes and showed that the r-mode evolution 
can progress in a number of different directions depending on unknown properties of the 
viscosity, leading to very complex consequent mode evolution and the associated GW signal. 

On the other hand, differential rotation, first suggested by Rezzolla et al. (2000, 2001a, 
2001b), is an unavoidable feature of nonlinear r-modes (Stergioulas & Font 2001; Lindblom 
et al. 2001; Sa 2004). Small values of a mentioned above are also supported by studies on 
the role of differential rotation, causing large scale drift of fluid elements, in the nonlinear 
evolution of r-modes (Sa & Tome 2005). In particular they parametrize the initial amount 
of differential rotation by K and then relate the parameter K to the largest amplitude that 
the r-mode can grow to. In this paper we will use the characteristic GW amplitude given by 
Sa & Tome (2006) to account for the average source spectrum. The adopted GW amplitude, 
parametrized by parameter K, in fact scales with the saturation amplitude a. Below we will 
discuss the influence of this quantity on the r-mode background. 

In addition to the average source spectra, the properties of AGWB also depend on GW 
source formation rate. Studies of cosmic star formation rate (CSFR) allow estimation of the 
birth rate of NSs. In the last decade our knowledge of cosmic star formation has been greatly 
improved due to advances in astronomical observation and hydrodynamic simulation. Here 
we take into account both observational and simulation-based CSFR models to obtain the 
NS formation rate and discuss their effects on the resultant GW background. In particular, 
we will investigate the role of the maximal redshift of different CSFR models in our results. 
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The high-frequency window of GW spectrum (10 Hz < f < a few kHz) is open today 
through pioneering efforts of the first-generation terrestrial IFOs, such as Laser Interferome- 
ter Gravitational Wave Observatory (LIGO) (Abramovici et al. 1992) in Livingston (LIGOL) 
and in Hanford (LIGOH), Virgo (Caron et al. 1997) near Pisa, GEO600 (Liick et al. 1997) 
in Hanover and TAMA300 (Ando et al. 2001) at Tokyo. Although GWs were not detected, 
an observational upper limit (^gw < 6.9 x 10"^) was placed on the energy density of SGWB 
at around 100 Hz, exceeding previo us indirect limits fro m the Big Bang Nucleosynthesis and 



the cosmic microwave background (lAbbott et al. 20091 ). In the future the SGWB from NS 



r-mode instability, among others, may offer an important detection target for the proposed 
second and third generation detectors represented by advanced LIGC10 (or advanced Virgcl^) 
and the Einstein Telescope (ET§|) respectively. 

Detectors throughout the world can act as a network in order to improve the detection 
ability to the SGWB. Two approaches of combining 2N IFOs are proposed in Allen & Romano 
(1999): (i) correlating the outputs of a pair of IFOs, then combining the multiple pairs, and 
(ii) directly combining the outputs of 2N IFOs. For any given rea£| IFOs it is necessary 
to compare these two optimal approaches of detecting the SGWB. Cella et al. (2007) has 
shown that the approach of combining multiple pairs of IFOs using Virgo, LIGO and GEO 
can improve the detection ability to the SGWB by simulating an isotropic GW background 
with an astrophysically-motivated spectral shape. Fan & Zhu (2008) compared the detection 
ability of the two approaches for stochastic GWs from string cosmology. 

Detectability of the r-mode background are demonstrated here by calculating signal-to- 
noise ratios (SNRs) for pairs of currently operating IFOs and advanced detectors at their 
design sensitivities. We also consider two approaches of combining 4 real IFOs to examine 
how many improvements can be obtained. The organization of this paper is as follows. In 
Section 2 we review works on the determination of CSFR and present five CSFR models. In 
Section 3 we derive the NS formation rate as a function of redshift using the adopted CSFR 
models. Then by combining the source formation rate from Section 3 and the characteristic 
GW amplitude of individual events, spectral properties of the r-mode stochastic background 



^http://www. ligo.caltech.edu/advLIGO/ 
^http://wwwcascina. virgo.infn.it/advirgo/ 
■^http://www.et-gw.eu/ 

^Here by "real" we mean considering the real overlap reduction functions for different detector pairs other 
than assuming an optimal value of unity since this function plays a crucial role in determining the frequency- 
dependent sensitivity of each detector pair to the stochastic background (Finn et al. 2009). Results in the 
present paper also show that using an optimal value of unity can lead to overestimates in signal-to-noise 
ratios of more than one order of magnitude. 
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are investigated in Section 4. We will discuss the detectability of r-mode background in 
Section 5 and finally Section 6 is devoted to our conclusions. 

Throughout the paper, the so-called 737 ACDM cosmology is assumed with Ho — lOO/i- 
kms-^Mpc~^ with h = 0.7 and = 0.3, Qa = 0.7 (e.g., Komatsu et al. 2009). 

2. Cosmic star formation rate 

The CSFR, which has tight connection with GW event rate, is of intense interest to 
many fields of astrophysics. For many years effort has gone into studying the cosmic star 
formation history (see, e.g., Madau et al. 1996; Hopkins 2004; Wilkins et al. 2008). Since 
CSFR is not a directly observable quantity, usually the rest-frame ultraviolet (UV) hght 
is considered to be an indicator (see Calzetti 2008 and references therein for details about 
CSFR indicators) of star formation because it is mainly radiated by short-lived massive 
stars. With the help of Hubble Space Telescope (HST) and other large telescopes, galaxy 
luminosity density of rest-frame UV radiation is studied, and then converted into CSFR 
density through the adoption of a universal stellar initial mass function (IMF) to calculate 
the conversion factor. Many authors have developed parameterized fits to the expected 
evolution of the CSFR with redshift. Firstly, following Porciani & Madau (2001), we adopt 
three different forms which model the CSFR density for redshifts up to z fa 4: 

p4z)i = 1.67Cihe5F{z)Gi{z) Meyi-'Upc-^ (1) 

with i=l, 2, 3 denoting the different models, a constant, /igs = h/0.65, Gi{z) a function 
of z and F{z) = [0^(1 + zf + f^A]^/V(l + ^f^'^- ^ constant factor of 1.67 is applied to 
account for conversion of a Salpeter IMF with a lower cutoff from 0.5Mq to O.IM0 (the one 
we will use below). F{z) convert the assumed cosmology from an Einstein-de Sitter universe 
to the ACDM cosmology (Porciani & Madau 2001). The first fit (hereafter SFRl) is given 
by Madau & Pozzetti (2000), with Ci = 0.3 and Gi{z) = e^-"^' / {e^-^^ + 45) where the CSFR 
increases rapidly from 2; = to reach a peak at around z = 1.5 and then gradually declines 
at higher redshifts. The second one (SFR2) is from Stcidcl et al. (1999) with C2 = 0.15 and 
G2{z) = e^'"^^ / {e^''^^ + 22) where the CSFR remains roughly constant at 2; > 2. The third 
model (SFR3) from Blain et al. (1999) has C3 = 0.2 and G3{z) = e3-05^-o-'7(e2 -«^ + 15) 
where CSFR increases at higher redshifts to account for effects of dust extinction. 

With the improvement in measurements of galaxy luminosity functions at a broad range 
of wavelengths, star formation history can be traced to higher redshifts. Here we consider 
the work by Hopkins & Beacom (2006), who refined the previous models up to redshift 
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z ~ 6 from new measurements of the galaxy luminosity function in the UV (SDSS, GALEX, 
C0MB017) and far-infrared (FIR) wavelengths (Spitzer Space Telescope). A parametric fit 
(hereafter HB06) is given by: 

. , , , 0.017 + 0.13Z 1., s 

P^(^) = S + (,/3.3)5.3 Meyr-^Mpc-^ (2) 



assuming 737 cosmology and a modified Salpeter A IMF (IBaldry fc Glazebrook 20031 ). Al- 
though the IMF used to derive HB06 is different from the standard Salpeter's, this will not 
introduce considerable errors to our results because the evolution of CCSNe rate based on 
the CSFR and on an assumed universal IMF is largely independent of the choice of the IMF 
(Madau 1998). 

Many authors have addressed the issues of calibrating the high-z CSFR through long- 
duration gamma-ray bursts (GRBs) (see, e.g., Yiiksel et al. 2008 and Kistler et al. 2009). 
Recently Wang & Dai (2009) uses latest GRBs data to constrain the CSFR up to 2 = 8.3. 
Meanwhile there are other methods to determine the high-z CSFR, such as observations of 
color-selected Lyman break galaxies (LBGs) (Bouwens et al. 2008) and Lya emitters (Ota 
et al. 2008). However, such calibrations cannot reach considerable agreements except for an 
overall decline at 2; > 4 (see Figure 1 of Yiiksel et al. 2008, Kistler et al. 2009 and Wang & 
Dai 2009). Due to huge uncertainties and the incompleteness of data sets we will not include 
them here. 

On the other hand, Springel & Hernquist (2003) derive the CSFR from hydrodynamic 
simulation of structure formation in ACDM cosmology. They study the history of cosmic 
star formation from the "dark ages", at redshift 2; = 20 to the present. The CSFR obtained 
in their study is broadly consistent with measurements given observational uncertainty and 
can be remarkably well- fitted by the following form (hereafter SII03): 

. P e-KY>[a{z - Zm)] 

PAz) = p.m- -, (3) 

p - a + a ■ exp[p[z - Zm)\ 

where a = 3/5, /3 = 14/15, Zm = 5.4 marks a break redshift and pm = 0.15 M0yr~"^Mpc~^ 
fixes the overall normalization. It is worth mentioning that they consider a ACDM model 
with the same parameters with our assumed 737 cosmology. 

In Fig. 1 we plot the CSFR predicted in the above five models as a function of redshift. 
SFRl, SFR2 and SFR3 show distinguishable features at z >2. SH03 peaks at a much higher 
redshift, between z = 5 and z = 6, than observation-based models (around z ~ 2). The 
cutoff of each curve in Fig. 1 corresponds to maximum redshifts of CSFR models: z^, = 4 
for SFRl-3, 2;^ = 6 for HB06 and = 20 for SH03. What Fig. 1 illustrates is our poor 
understanding about star formation history at high redshifts from astronomical observations. 
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Note that there are some other CSFR models, similar to or different from the five 
models adopt here, not included since our aim is not to make a complete survey on this 
issue but to phenomenologically investigate its influences on the SGWB from an ensemble of 
astrophysical sources. We refer readers to Calura & Matteucci (2003), Daigne et al. (2004), 
Bromm & Loeb (2006), Nagamine et al. (2006) and Fardal et al. (2007) for details of other 
studies on the determination of CSFR. In the following sections we will investigate how 
different CSFRs affect the rate of NS formation and spectral properties of AGWB. 



Since the evolving rate of CCSNe closely tracks the star formation rate, using CSFR 
models presented in Section 2 we can estimate the number of NSs formed per unit time 
within the comoving volume out to redshift z (Ferrari et al. 1999a): 



where /?*(;z) is the CSFR density, dV/dz is the comoving volume element, and $(m) is the 
IMF. Here we assume that each CCSN results in either a NS or a BH and take a NS progenitor 
mass range of 8Mq — 25M0. In order to make comparison with Ferrari et al. (1999b) we 
also consider a lower upper limit for NS progenitor masses mmax = 20Mq as indicated from 
core collapse simulations by Fryer (1999). However according to Belczynski & Taam (2008) 
the mass of NS progenitor might be greater than 4QMq for stars in a binary system. So we 
will also include a higher limit of mmax = AQMq for our calculations of NS formation rate 
(In Sigl 2006 the progenitor mass to form a NS ranges from lOM© to AQMq). 

Note that in some studies (e.g.. Coward et al. 2001; de Araujo et al. 2004; Regimbau & 
Mandic 2008) with respect to AGWB there is an additional (1 + z) term in Eq.(|l]) dividing 
the CSFR to account for the time dilatation of CSFR due to cosmic expansion. Here we 
do not include such a term according to de Araujo & Miranda (2005) who argue that the 
inclusion of this additional term is inadequate. 

To integrate through Eq.(|4]) one still needs to know the forms of dV/dz and <I>(m). 
Following Regimbau & Mandic (2008), the comoving volume element is related to z through 



3. 



Neutron star formation rate 




dz 




(5) 



Ho E{n, z) ' 



where Hq the Hubble constant, E{Vl, z) = a/^a + f2m(l + zY and r[z) the comoving distance 
related to the luminosity distance by di = r^i)- + z). 
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We consider the standard Salpeter IMF: $(m) = Am^^-^'^^'> with x = 1.35, where A is a 
normahzation constant, obtained through the relation m^{m)dm = 1 with mi = O.IMq 
and m„ = 125Mq. Then we plot the NS formation rate R^si^) defined in Eq.(jl]) for the 
CSFR models presented in Section 2 with a modest value of mmax = 25M0 in Fig. 2. Note 
that SFRl, SFR2 and SFR3 show quite different behaviors for z > 2.5 and observation-based 
models give rise to more NS formation than SH03 up to respective redshift limits, but SH03 
predicts a much higher cumulative NS formation rate for z > 10. 

In Table 1 we present the total number (per unit time) of CCSN explosions leaving 
behind a NS out to corresponding redshift limits for the five CSFR models, and for three 
values of mmax : 20Mq, 25Mq and AOMq. We compare the results obtained here with those 
in Ferrari et al. (1999b), and find a factor of ~ 2 — 3 enhancement for the total NS formation 
rate, which is mainly due to differences of CSFR models and cosmology terms (e.g., different 
forms for the comoving volume element). 



4. Spectral properties of the SGWB from NS r-mode instability 

In this section we will evaluate the spectral properties of the stochastic background 
produced by an ensemble of newly born NSs with nonlinear r-mode instabilities. Initially, 
let us review the formalism used to characterize the AGWB. 

It is useful to characterize the spectral properties of a SGWB by specifying how the en- 
ergy is distributed in frequency domain. Explicitly, one introduces a dimensionless quantity, 
f^GW given by: 

^Gw('^obs) = — , (6) 

pc d In z/obs 

where pcw is the GW energy density, u^hs the frequency in the observer frame and pc = 



Table 1: Total rate of NS formation in five CSFR models for three values of the upper limits 
of NS progenitor masses. 



Model (redshift limit) 


?^max = 2OM0 


25Mq 


4OM0 


SFRl (2^ = 4) 


30.0 


33.2 


37.4 


SFR2 (2^ = 4) 


39.3 


44.5 


49.1 


SFR3 (-2^ = 4) 


47.2 


52.2 


58.9 


HB06 {z^ = 6) 


47.5 


52.6 


59.3 


SH03 {z^ = 20) 


62.0 


68.6 


77.4 
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3Hq/87cG is the critical energy density required to close the universe today. For a stochastic 
background of astrophysical origin, the energy density is given by: 

^Gw(i^obs) = ^r^-Pl(^'obs), (7) 

where the spectral density of the flux at the observed frequency Uobs is defined as 

Fu = J fu{i^ohs)dR, (8) 

where /^(i^obs) is the energy flux per unit frequency (in erg ■ cm~^ ■ Hz^^) produced by a 
single source and dR is the differential GW event rate. 

The energy flux per unit frequency fu{i^ohs) can be written as follows (Carr 1980) 

/.(^obs) = ^hl, (9) 

where he is the dimensionless amplitude produced by an event that generates a signal with 
observed frequency Uohs- 

In order to obtain the spectral properties (e.g., the values of flaw as a function of Uohs) of 
the r-mode stochastic background, we have the differential rate of source formation dR^s{z) 
through Eq.(j4l) and still need to know the energy flux emitted by a single source. 

It has been shown that differential rotation can significantly influence the detectability 
of GWs emitted by a spinning-down newborn NS due to r-mode instability (Sa & Tome 
2006). Studies of mode-mode coupling in rotating stars also indicate that the maximum 
amplitude that r-mode can grow to is much smaller than previously estimated (Arras et al. 
2003). Here we use the characteristic GW amplitude given by Sa & Tome (2006): 

, , , 5.5 X 10-22 ri^ f20Mpc\ , , 

where i^' is a constant giving the initial amount of differential rotation associated with the 
r-mode, and lies in the interval —5/4 < K < 10^'^ (see Sa & Tome (2006) for details), 
= ^ohs{^ + z) is the frequency in the source frame, i^max the maximum frequency of emitted 
GWs given by 2f2K/37r where Qk = 5612 Hz is the Keplerian frequency at which the star 
starts shedding mass at the equator, assumed to be the initial value of the angular velocity of 
the star, and d^ is the luminosity distance to the source. Note that the saturation amplitude 
a oc {K + 2) -1/2 (Sa & Tome 2006), which means the GW amplitude is proportional to a, 
same as those given in Bondarescu et al. (2009). In the following calculations we should keep 
in mind that the parameter K used in this paper is equivalent to the saturation amplitude 
a of NS r-mode instability. 
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From the above equations we can obtain the dimensionless energy density: 



O _ 4vr^(l.l X 10-^°)^ ^ 



/ / p4z){1 + z)[—^] —^{m)dmdz 
Jz^^„ Jm^;„ \ di ^ dz 
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Thus, by setting a value for K one can calculate ^gw numerically through Eg. ( ITT]) 
combined with corresponding equations for CSFR, comoving volume element and IMF. Here 
we set minin = 8M0 and mmax = 25M0, while ^min and Zmax can be determined in such a 
way: since frequencies of emitted GWs in the source frame range from z/min = 77 — 80 Hz to 
^'max = 2nK/37r = 1191 Hz, where the minimum frequency corresponds to the final angular 
velocity of the star - 0.065f2K for K = —5/4 and 0.067f2K if -ft" 3> 1, we have Umin/ + z) < 
^ohs < ^lasLx/ (1 + ^)! which means sources with different redshifts that produce a signal at the 
same frequency i^obs should meet the condition: z/min/i/obs — 1 < -2 < ^'max/^'obs — 1- Besides, 
we consider signals emitted at early epochs up to the present {z > 0) and take into account 
the maximal redshift (z^,) of CSFR model. Then we obtain Zmin = iiiax(0, z/min/z/obs — 1), 
^^max = min(z*, z/jnax/^'obs ~ 1)) which is similar to that of Owen et al. (1998) where 2;* ~ 4 is 
considered to be the maximum redshift where there was significant star formation. 

In Fig. 3 we plot the dimensionless energy density ^gw calculated for the five CSFR 
models presented in Section 2 by setting K at its minimal value: K = —5/4 corresponding to 
the smallest amount of differential rotation at the time when the r-mode instability becomes 
active. However, as emphasized by Sa & Tome (2005), if K is small, namely, 0, it is 
necessary to consider other nonlinear effects like mode-mode couplings in the calculation of 
a, which will again limit the maximum r-mode amplitude to values much smaller than unity 



(lArras et al. 20031 ). In this respect, our choice of a minimum K results in an unrealistically 



high upper limit for r-mode background. 

It is worth noting from Fig. 3 that no obvious differences are recorded for the three 
curves of SFRl, SFR2 and SFR3, and observation-based CSFR models give rise to stochastic 
backgrounds about two times stronger than that of SH03 over a broad frequency band, 
although SH03 leads to a much higher NS formation rate. The sharp contrast between Fig. 
3 and Fig. 2 indicates that the main contribution to the GW background comes from low- 
redshift sources because those events happened at higher redshifts have minor influences due 
to the inverse squared luminosity distance dependence of the single event energy flux. For 
the same reason, poor observational understanding of high-z star formation history (see Fig. 
1) is not severe to studies about AGWB here. 

To assess the role of in our results, we choose SFR2 (since it remains constant for 
z > 2), set three values for z^, {z^ = 4,10 and 20), and then plot the ^qw in Fig- 4. It 
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is surprising that the three curves exhibit almost the same pattern in the frequency range 
^'obs > 100 Hz, and extending the redshift hmit from 4 to 20 results in a growth of the 
lower-frequency background. Increasing 2;* will enhance the background at lower frequencies 
and even enable some formerly "unavailable" low-frequency signals to emerge. This low- 
frequency GW "tail" can be accounted for by the contribution from high-redshift sources. 
However if CSFR is much lower at high-z, this effect will be negligible. Thus Fig. 4 further 
support the conclusions from Fig. 3 and indicate that the most significant contribution to 
an AGWB comes from GW events occurring at redshifts z < 4. 

From Eq. f fTTj) we find that ^gw depends on the values of z/max ~ and K {flew oc 
~ a^). As suggested by Ferrari et al. (1999b), the Keplerian velocity Qk may have 
a broad distribution due to different masses and radii of rotating NSs. Here we arbitrarily 
set z/jnax ranging from 1000 Hz to 2000 Hz. For K we set —5/4, 100 and 10^ corresponding 
to a = 1,0.1 and 10^^ respectively. Then we adopt HB06 as the CSFR model (below we 
will use only HB06) and plot f^cw as a function of observed frequency for different values 
of a and i^max in Fig. 5. We can see a higher peak for Qqw when we increase the maximum 
emitting frequency, while the r-mode background for smaller z/max is slightly enhanced at 
lower frequencies (< 400 Hz). On the other hand, increasing the amount of differential 
rotation significantly reduce the closure density and then affect the detectability of r-mode 
background as we will discuss later. Considering a maximum value of 10~^ for a, a realistic 
estimate of r-mode background should have a energy density at most ~ 10^^^ like the lowest 
curve shown in Fig. 5. 

Another important quantity of the AGWB is the so-called duty cycle, which classifies 
the stochastic backgrounds in terms of continuous background, popcorn noise and short noise 
(Coward & Regimbau 2006): 



where f is the average time duration of the GW emission from a single source at the source 
frame, which dilated to f{l + z) by the cosmic expansion, and (ii?Ns(^) is the differential rate 
of NS formation in Eq.(jl]). It has been suggested that differential rotation can remarkably 
influence the long-term spin and thermal evolution of NSs by prolonging the duration of the 
r-modes (Yu et al. 2009). In view of the prolonged r-mode, the spinning-down phase can 
last even longer than 1 yr, which indicates a duty cycle >~ 10^. In the next section, we will 
discuss the detectability of this continuous GW background. 




(12) 
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5. Detectability 



5.1. Detecting the r-mode background with a network of IFOs 



We cannot reach sufficient sensitivity for detection of a SGWB with a single terrestrial 
IFO since the output of a detector is dominated by the noise rather than by the signal due 
to the stochastic background itself (Allen 1996a; Maggiore 2000). The optimal strategy to 
search for a SGWB is to cross correlate measurements of two or more detectors. It has 
been shown that after correlating signals of two detectors for a (i) isotropic, (ii) unpolarized, 
(iii) stationary, and (iv) Gaussian stochastic background (see Allen & Romano 1999 for 
discussions of these assumptions), the optimal SNR during an integration time T (here we 
assume T = Syr ~ 10^ s) is given by an integral over frequency /: 



where -Pi(/) and -P2(/) are the power spectral noise densities of the two detectors and 7(/) 
is the so-called overlap reduction function, ffist calculated by Flanagan (1993). This is a 
dimensionless function of frequency and determined by the relative locations and orientations 
of two detectors. For 7(/), we refer readers to Flanagan (1993), Allen & Romano (1999), 
Maggiore (2000) for more details. 

To assess the detectability of the r-mode background, we will calculate the SNRs for 
several pairs of detectors for ficw computed with HB06 CSFR model, K = —5/4 and 
^max = 1191 Hz (unless otherwise stated we use these parameters in Section 5). Here we 
consider the four IFOs which are in routine operations - LIGOH (4km), LIGOL (4km), Virgo 
(3km) and GEO (600m), as well as the second generation detectors - advanced LIGO and 
advanced Virgo. Design sensitivity curves of these detectors are shown in Fig. 6. It is worth 
mentioning that the ffist-generation GW interferometric detectors have taken data at, or 
close to, their design sensitivities (Fairhurst et al 2009). In particular, the design sensitivity 
curve for initial LIGO was almost attained by its S5 run. In the following calculations we 
will use real 7(/) for different pair^ of IFOs unless otherwise stated. 

Two approaches of combining 2N detectors to improve the detection ability to the SGWB 
are proposed in Allen & Romano (1999): (i) correlating the outputs of a pair of detectors, 
then combining multiple pairs (combining pairs, "c-p"), and (ii) directly combining the 
outputs of 2N detectors (directly combining, "d-c"). For the ffist approach, the squared 




(13) 



^Data of locations and orientations of are taken from Allen (1996b). 
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SNR is given by: 



and for the second one: 



: 



^ ^ (12) f ^\ (34) f ^\ . . . (2Ar-l,2Ar) f 



Nlopm ViVy \N} \N 



+ all possible permutations . (15) 



We show in Table 2 the SNRs calculated for different pairs of LIGOH, LIGOL, Virgo 
and GEO, and for two combinations of these four IFOs. We consider here two cases rep- 
resenting two real networks of first/second generation IFOs: Case 1 is for these four IFOs 
with design sensitivities; Case 2 consists of two advanced LIGO detectors and two advanced 
Virgo detectors both with proposed sensitivities. Note that SNRs in Table 2 are lower than 
unity even for pairs of advanced detectors. The most promising one {SNR = 0.58) comes 
from combining pairs of four advanced IFOs. If we assu me an optimized value o f unity for 



7(/), which is only possible for co-located GW detectors ( iFotopoulos fc LSC 20081 ). the SNR 
is 11.0 and 4.1 for a pair of advanced LIGO and advanced Virgo IFOs respectively. We note 
that by considering new detectors with comparable sensitivities to advanced LIGO, such as 
LCGT in Japan (Kuroda et al. 1999) and AIGO in Australia (Blair et al. 2008), it could be 
possible to reach a higher, but still not significant SNR with a network of second-generation 
IFOs. 

In order to obtain some detectable parameter space we need at least one order of mag- 
nitude higher SNRs than those in Case 2 of Table 2. Then we reduce the noise power 
spectral densities of advanced LIGO and advanced Virgo by a factor of 10 by hand, and 
investigate the role of differential rotation (K) and maximum emitting frequency (z/max) in 
the detectability of the r-mode background. We are motivated here by the fact that third- 
generation detectors like ET can reach a sensitivity roughly an order of magnitude better 
than that of advanced LIGO (Hild et al. 2008). 

In Fig. 7 we plot the SNR as a function of K for H-L, H-V and L-V pairs. As a natural 
result from Fig. 5, the detectability of SGWB from r-mode instability is drastically reduced 
to as if approaching 10. The higher SNR of H-L pair reflects the lower noise level of 
advanced LIGO. Due to similarity of the overlap reduction functions (see Fig. 2 of Fan & 
Zhu 2008) no significant difference is shown between L-V and H-V pairs. 

Fig. 8 shows the SNR evolution with z/max for H-L, H-V and L-V pairs. In contrast 
to Fig. 5, in which a higher peak value of ^gw was obtained for larger z/max, there are no 
identical features for SNR evolutions here. For H-L pair, we note that SNR varies inversely 
as the increase of I'rna.x- This can be explained that the low- frequency GW background is 
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enhanced for smaller i/max while the growth of high-frequency background due to larger t'max 
is suppressed by the 1//^ term in Eq. f lT3|) . On the other hand, we can reach higher SNRs 
for H-V and L-V pairs for larger z/max (> 1600 Hz). This unique feature can be attributed to 
particular evolving behaviors of 7(/) for the two pairs (see again Fig. 2 of Fan & Zhu 2008 
and we will find that 7(/) is extremely close to zero at high frequencies for H-L, while it 
still fluctuates above and below zero up to 1000 Hz for H-V and L-V pairs). This is strongly 
supported by three curves in Fig. 9, where we plot the SNR as a function of i/max by assuming 
7(/) = 1 , show exactly the same evolving pattern as the curve of H-L in Fig. 8. 

For a detection rate 90% and a false alarm rate 10%, the total optimal SNR threshold 
should be 2.56. In order to evaluate the promise of detecting the r-mode background we 
present in Fig. 10 the regions in the (t'max; K) plane where SNR could be higher than 
2.56 for multiple detector pairs. It is shown that the detectable parameter space is quite 
limited even for a real network of third-generation IFOs. In particular, a strong constraint 
for {K < or equivalently a ~ 1) is obtained. Meanwhile we find that two approaches 
of combining multiple IFOs can effectively improve the detection ability as compared to 
H-L pair, and "c-p" method performs better than "d-c" approach in terms of detecting the 
r-mode background. 



5.2. The detection prospects of future detectors 

In this section we discuss further the detection prospects of a SGWB from NS r-mode 
instability by networks of third generation instruments. 

Fig. 10 indicates that only when the initial amount of differential rotation is near its 
minimum value {K ^ 0), the r-mode background could be detectable by a real network of 
third-generation GW detectors. However, as emphasized by Sa & Tome (2005), if K is near 
its minimum value, other nonlinear effects like mode-mode couplings should be included in 
the calculation of satura tion amplitude a . Still in this case, a will be limited to values 



much smaller than unity (lArras et al. 20031 ) . In this respect the physically reasonable values 



of parameter K introduced in Sa & Tome (2005) should be larger than 10^ i n order to be 



consistent with a maximum saturation amplitude a < 10 ^ (ILin fc Suen 20061 ). 



It seems most probable that the r-mode background is not going to be detectable. Given 
the relation SNR ~ f^cw ~ ^ 10~^ saturation amplitude (compared with a value of order 



®In fact such an improvement is negligible considering that SNR is extremely sensitive to the saturation 
amplitude of r-mode (oc a^). 
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unity) will reduce the SNR by 4 orders of magnitude for any measurements of SGWB from NS 
r-mode instability. In particular, to reach the same SNR requires 10000 times improvements 
in sensitivities of two detectors. 

If we leave the small a issue (the dominant aspect) aside, it is still quite pessimistic to 
detect the r-mode background as seen from Table 2 and Fig. 10. This is due to the following 
facts: 1) the closure density of this background peaks at much higher frequencies than the 
most sensitive frequency band of ground-based detectors; 2) the minimum frequency of r- 
mode GW signal z/min = 80 Hz just misses the frequency band (1-60 Hz) where the overlap 
reduction functions of real detector pairs are significant; 3) we should not forget that we 
assumed all NSs are born with angular velocities near their maximal value f^K. This is not 
necessarily always the case since it would make more sense to consider some fraction of NSs 
are born with rapid spins (Owen et al. 199 8). Actually it has been suggested that most 



NSs are born with very small rotation rates (jSpruit fc Phinney 19981 ). Current population 



synthesis stu dies favor spin periods of NSs at b irth in the range from tens to hundreds of 



milhseconds (lOtt et al. 2006 I : IPerna et al. 2008f l. 



Realistic overlap reduction functions are adopted here for multiple ground-based IFOs. 
Recall that two co-located advanced LIGO detectors (7(/) = 1, SNR=11) perform even 
better than combining four IFOs with 10-fold better sensitivities (SNR=5.8). With this 
in mind we also evaluate the detectability of r-mode background with ET, assuming two 
detectors located in Cascina, of triangul ar shape (60° between the two arms) and separated 



by an angle of 120° (iHowell et al. 20101 ). The 7(/) benefits a lot from this configuration, 
nearly linearly decreasing from —0.372 at 1000 Hz to —0.375 at 1 Hz. We adopt the ET-B 
sensitivity from Hild et al. (2008). For 3- year integration, a SNR of 2.56 requires K <~ 150 
corresponding to a ~ 0.1. 

Then we convert the constrains on i^T or a to absolute numbers: the total emitted GW 
energy associated with NS r-mode instabihty that enable the stochastic background to be 
detectable by future detectors. The energy flux of single event can also be written as: 

where dEcw / dv is the gravitational spectral energy. Combined with Eq. ([9]) and Eq. ffTOj) we 
can obtain the spectral energy density and thus the total GW energy emitted by individual 
sources that contribute to this background. We give the results (required energy level in 
order to obtain a SNR of 2.56 by 3-year integration) as follows (in MqC^): i) 1.8 x 10~^ for 
a real network of "third-generation" IFOs; ii) 9.4 x 10~^ for two co-located advanced LIGO 
detectors; iii) 2 x 10^^ for two detectors with ET-B sensitivity. 
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6. Conclusions 

We revisit the possibility and detectability of a SGWB produced by a cosmological 
population of young NSs with an r-mode instabilities. Source formation rate is accounted 
for by using a set of CSFR models, both observational and simulated. Our results show that 
the resultant GW background is insensitive to the choice of CSFR models (although they 
predict quite different NS formation rates) , but dependent on the evolving behavior of CSFR 
at low redshifts {z < 2). This is in good agreement with that in Howell et al. (2004). But 
here we further investigate the effect of the maximal redshift of CSFR models on the AGWB 
and find that high-redshift {z > 4) sources could form a GW "tail" in the low-frequency 
side (< 30 Hz, this number depends on particular source spectrum and thus only applies 
to r-mode background) and has no effect on the high-frequency background. Such an effect 
will be negligible if high-z CSFR is much smaller. Our Fig. 3 and Fig. 4 indicate that the 
most significant contribution to the GW background of astrophysical origin comes from GW 
events occurring at redshifts z < 4. 

The characteristic GW amplitude parametrized with the initial amount of differential 
rotation (K) during r-mode evolution is adopted as the average GW signal for individual 
sources. While a minimum K corresponding to a saturation mode amplitude a ~ 1 is 
assumed, the energy density JIgw has a maximum amplitude at around 3 x 10~^, agreeing 
well with those of Owen et al. (1998) and Ferrari et al. (1999b). However since we know 
that the maximum amplitude a that r-mode can grow to is at most lO^'' — 10~^. This means 
the physically reasonable values of parameter K is at least 10'^. Consequently a realistic 
estimate of i^ow for r-mode background should be at most ~ 10~^^. 

We further consider multiple IFOs, including the first-generation ones that are in routine 

operations and upgraded counterparts, for the detection of r-mode background. We also 
illustrate how a network of ground-based IFOs could improve the detection ability to this 
GW background. Since the detectability is dominated by the square of saturation amplitude 
{SNR oc a^), it is likely that the r-mode background will not be detectable for any future 
detectors. Still we give the constraints on the total emitted GW energy associated with 
this mechanism to enable a detection with 90% detection rate {SNR — 2.56) by 3-year 
cross correlation as: ~ 10~^ MqC^ for two co-located advanced LIGO detectors and 2 x 10~^ 
for two IFOs with ET sensitivity. Considering the relatively certain NS formation rate, 
these constraints might be applicable to alternative emission mechanisms associated with 
NS oscillations and instabilities. This requires further investigation. The requirement on 
GW energy level could be lower if more signals are emitted near the most sensitive frequency 
band of ground-based detectors (~ 40 — 200 Hz). Through reasonable assumptions of average 
source spectra, the lowest detectable (in terms of stochastic background) GW energy for 
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individual source will be 10~^ MqC^ for third generation detectors like ET ( iZhu et al. 2010l ). 
Overall, for the detection of SGWB from NS instabilities, more efficient emitters are required 
(see, e.g., Andersson et al. 2010 for reviews of GW emission from NSs and Kastaun et al. 
2010 for details of a possible more efficient mechanism - f-mode). 

While the SGWB from NS r-mode instabili ty is difficult to detect, the associated GW 
signal is still detectable in terms of single events ( lOwen 2010 I ), although this possibihty also 
depends strongly on the saturation amplitude. For instance, initially it was believed that 
GWs from r-mode instab ility in a newborn NS co uld be detected by advanced LIGO out 
to a distance of 20 Mpc (lOwen &: Lindblom 20021 ). However even for the most optimistic 
case in Bondarescu et al. (2009), the detectable distance for advanced LIGO is only 1 Mpc. 
The LIGO Scientific Collaboration and Virgo Collaboration have already performed many 
searches for periodic GWs from rapidly rotating NSs inclu ding the first search targeting 
the youngest known NS - Cassiopeia A (lAbadie et al. 20101 ). Although no GWs have been 
detected, direct upper limits on GW emissi on from known puls ars like the Crab pulsar have 
beaten down the indirect spin-down limits ( lAbbott et al. 20081 ). 
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Table 2: Concerning the detection of the r-mode background, we present the SNRs calculated 
for multiple detector pairs of LIGOH (H), LIGOL (L), Virgo (V) and GEO (G) and two 
approaches of combining the four IFOs - combining pairs (c-p) and directly combining (d-c). 
Case 1 corresponds to first-generation detectors and Case 2 adopts sensitivities of advanced 
LIGO and advanced Virgo, while G refers to an assumed interferometer with the same 
sensitivity as advanced Virgo at the GEO site. 



Case 


H-L 


L-V 


H-V 


V-G 
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1.8 X 10^3 


1.0 X 10-3 


8.5 X 10-^ 


2.6 X 10-^ 
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0.40 
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0.17 


0.25 


Case 


L-G 


H-G 


c-p 


d-c 
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1.4 X 10"^ 


9.7 X 10-5 


2.3 X 10-3 


5.0 X 10-^ 
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0.16 


0.11 
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Fig. 1. — Evolution of the CSFR density p^{z) (in MqJ"^ ^Mpc predicted in five param- 
eterized models (see text). 
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Fig. 2. — Evolution of the NS formation rate under different CSFR models. 
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Fig. 3. — The dimensionless energy density ficw as a function of observed frequency fobs, 
calculated for five CSFR models and by setting = — 5/4, i/max = 1191 Hz. 
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Fig. 4. — ^GW as a function of observed frequency z^obs, calculated for SFR2 model with 
three different values of maximum redshift z^, and for other parameters same as Fig. 3. 
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Fig. 6. — The designed noise power spectrum of the GEO, Virgo, initial LIGO, advanced 
Virgo and advanced LIGO. Data are taken from Lazzarini et ah 1996, Punturo 2004, 
http : / / www, aei . mpg . de /j rsmith /geo curves . html | |http : / / wwwcascina . virgo . inf n . it / advirgo / 
and https:/ /dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=2974. 
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Fig. 7. — SNR as a function of K for H-L, H-V and L-V pairs when we set Umax = 1191 Hz. 
We assume a factor of 10 folds improvement in sensitivity for advanced LIGO and advanced 
Virgo detectors. 
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Fig. 9. — SNR as a function of z/max assuming 7(/) = 1 for two advanced LIGO (adL) 
detectors, two advanced Virgo (adV) detectors and the combination between them (here no 
improvement of sensitivity is assumed). 
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Fig. 10. — The region of parameter space for which SGWB produced from an ensemble of NSs 
with r-mode instabihty is detectable by H-L pair (here an order of magnitude improvement for 
advanced LIGO sensitivity is assumed), the approach of directly combining (labeled "d-c") 
four third-generation detectors (an order of magnitude improvement in detector sensitivities 
for Case 2 in Table 2 is assumed) and the approach of combining multiple pairs of detectors 
(labeled "c-p"). With 10% false alarm and 90% detection rate, the region above the curves 
shows the detectable parameter space. 



